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Abstract 

The problem of the Hamiltonian matrix in the oscillator and or- 
• thogonalized Laguerre basis construction from a given S-matrix is 

treated in the context of the algebraic analogue of the Marchenko 

(N 
O 

Q 

00 

The J-matrix pQ theory of scattering is based on the fact that the £th partial 
wave kinetic energy or the Coulomb Hamiltonian H° is represented in a cer- 
tain square-integrable basis set by an infinite symmetric tridiagonal matrix. 

\D ' 
O 

Q\ ' problem for H° can be solved analytically. The J-matrix method yields an 

O ' 



method. 



1 Introduction 



exact solution to a model scattering Hamiltonian where the given short-range 

(- . N-l 

potential is approximated by truncating in a finite subset yF n j ■ 

In Refs. [2JIH] an inverse scattering formalism within the J-matrix method 
has been proposed, where the matrix ||V„ im || of the potential 

c3 ■ N-l 

Vi(r, r') = hw J2 tiS?) Vn,m<Pl(x') ( 1 ) 



n, m=0 



with the oscillator form factors 

= (-i) 



2n\ 



x e 



1.-^/2^+1/2^2) (2) 



is determined from a given S-matrix. Here, x = r/p is the relative coordinate 



in units of the oscillator radius p = yh/puj, p is the reduced mass. 

Obviously a correlation can be made between the J-matrix method and 
a discrete model of quantum mechanics, within of which a finite-difference 
Schrodinger equation is used. As a result the J-matrix versions of the 



^This work has been done partially while the author was visiting the Institute for 
Nuclear Theory, University of Washington. 
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Gel'fand-Levitan-Marchenko method algebraic analogue can be formulated. 
For instance, the J-matrix method formally and computationally is quite 
similar to the R-matrix theory. It is this analogy that the previous J-matrix 
version of inverse scattering theory [2JIHJII] [also see [H]] leans upon. Within 
the J-matrix approach the discrete representation of the Green function in 

finite subspace of the basis functions \ 4> n \ 



G( 



(el-h)- 



(3) 



is used. Here, I is identity matrix and h is the truncated Hamiltonian matrix 
of order N in the oscillator basis (j2J). We measure the energy E in the 
units of the oscillator basis parameter hu, i.e. E = hue and e = q 2 /2, 
where q is the dimensionless momentum: q = hp. In particular, the element 



[G(e)] 



N-1,N-1 



Vn(c) can be presented in the two rational forms [21 E]: 



JV-2 

n (e - Hj) 

= 15 

n (e - \j) 

j=0 



(4) 



where {\j}^ = Q and {^j}^ = Q satisfy the interlacing property 
\ < fi < \i < . . . < Xn-2 < < AiV-l, 



and [7 



V N {e) = E 



JV-l 72 

Z AT-l,j 

e-\j' 



(5) 



3=0 



Here, {Zn-ij} = q are the elements of the iVth row of the eigenvector matrix 



Z of the truncated Hamiltonian matrix h. The sets {Aj}^ 1 , {pj}j =u 

{Xj}-~y, {Zn-ij^Sq are derived from the S-matrix which is intimately 
connected with Vn [see Eq. (j3HJ)]. Notice that the both sets of the spec- 
tral parameters determine unique [apart from the off diagonal elements sign] 
Hamiltonian matrix h of a Jacobi form [HI Ej 



N-2 



and 



/ a 61 
bi ai 





b 2 

h 







63 

X 
X 



X 
X 

6jv— 1 



b N -i 
ajv-i / 



(6) 
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Hence, the sought-for potential matrix ||V^ m || is also of a Jacobi form. Recall 
that the kinetic energy operator 



matrix representation \\T* m \\ = £-\\H® m \\ in the harmonic oscillator basis 
(J2J) is symmetric tridiagonal [TJ [7] : 



Tl n = U2n + £ + l), 



Tl n+1 = Ti +1>n = -|^(n + l) ( n + < + |). 

Thus the inverse scattering problem within the J-matrix approach admits 
of the solution in the tridiagonal Hamiltonian matrix form. In this regard 
the J-matrix method is similar to a discrete model of quantum mechanics, 
in the framework of which the Hamiltonian matrix representation is also 
symmetric tridiagonal. Note that the tridiagonal matrix representation of 
both the kinetic energy operator and the Hamiltonian is fundamental for a 
finite-difference analogue of the Gel'fand-Levitan equations [see e.g. [TO]] as 
well as for a discrete version X 1 j of the Marchenko equations. As shown 
below, the similarity between the J-matrix method and a finite-difference 
approach can be plainly extended to the inverse scattering formalism as well. 
In the present paper, in particular, an inverse scattering J-matrix approach 
via the Marchenko equations (JME) is given. 

In JME the expansion coefficients c n of the wave function if) in terms 
of the L 2 basis set play a role similar to that of the values if> n = if>(x n ) at 
points x n = nA within the finite-difference inverse scattering approach. Here, 
the completeness relation for the solutions c n of the Schrodinger equation 
discrete analogue is also exploited which involves an integration c„c^ over 
the energy from zero to infinity. This raises the question as to whether the 
integrals converge. As shown below, taking account of the phase shift Sg(k) 
corresponding to the potential (0) of finite rank iV asymptotic behavior at 
large k provides the convergence of the integrals. Generally, with a potential 
((H) matrix of finite order iV it is possible to reproduce the phase shift Si 
only on a finite energy interval [0, eo] with eo < Ajv-i- This is why the 
eigenvalue Xn-i and the corresponding eigenvector component Zn-i,n-i are 
the variational parameters within the previous J-matrix version |3 13 El El 
of the inverse scattering theory. By contrast, in JME the phase shift, even 
if modified in accordance with the Sf asymptotic feature, on infinite energy 
interval is used. As a result JME has not variational parameters (apart from 
N and p). 
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The elements of the J-matrix method formalism are presented in Sec. 2. 
In Sec. 3 the inverse scattering J-matrix approach in the context of the 
Marchenko equations is formulated. The features of JME numerical real- 
ization are discussed in Sec. 4. In Sec. 5 JME is expanded to the Laguerre 
basis case. Here, we are dealing with the tridiagonal Hamiltonian matrix con- 
struction in an orthogonalized Laguerre basis set, in which the kinetic energy 
operator matrix is also tridiagonal. In Sec. 6 we summarize our conclusions. 

2 The direct problem 

The oscillator-basis J-matrix formalism is discussed in detail elsewhere. We 
present here only some relations needed for understanding the inverse scat- 
tering J-matrix approach. Within the J-matrix method, the radial wave 
function ip(k, r) is expanded in an oscillator function (J2J) series 

oo 

^(*,r) = £ ( 9 ) 

n=0 

In the assumption that the Hamiltonian matrix is of the form © the func- 
tions c n are the solutions to the set of equations 



a c (k) + &i d(k) =ec (k) 
b n Cn-i(k) + a n c n (k) + b n+1 c n+1 {k) = ec n (k), n = 1, 2, . . . . 
The asymptotic behavior of c n (k) for k > as n — > oo is given by 

c n (k) = f n (k) = ~ [ C { J(q) - S(k) C { n %q) 
Here, the functions 

= Cn,e{q) ± iS nJ {q), 

<Sn. 



(10) 



M-^r^^ qV2L ^% (12) 



CnAq) = ^f^Tf) e- q2/2 F(-n-£- 1/2, —t + 1/2; g 2 ) 

obey the "free" equations 

T n, n -i e n~i(q) +Tl n e n (q) + Tl n+1 e n+1 (q) = ed n (q), n = 1, 2, . . . . (13) 
S n> i satisfy in addition the equation 
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Besides, S n j meet the completeness relation 

oo 

2 r 

- / dkS ni e{q)S m . e {q) = <f n TO . (15) 

71 J 



Notice that 

~ oo 

S(r)= E S n>i (q) cf>i(x) 



n=Q 
oo 



(16) 



C(r)= EC M (# n (i), 

n=0 

subject to the asymptotic condition [Zj 

S(r) ~ sin(/cr — £tt/2), 

r— >oo 

C7(r) ~ cos(A;r -in/2). 

r— >oo 

As for the coefficients of the expansion 

oo 

Vv(0 = E^M^) (17) 

n=0 

of the normalized bound state wave function ip v with the energy — k„/2, 

Cn(iKv) = fnfav) = M v i* C { +}{\K v p) (18) 

holds as n — > oo. Here, .Mt, is the bound state normalization constant which 
is related to the residue of the S-matrix 



i Res S(k) = (-if M 2 V . (19) 

It can be easy verified that from the completeness relation for the solutions 

i[>(k, r), ip u (r) 



oo 

2 



- / dk ^(fc, x) i>{k, y) + 5>„0*0 = - 2/) (20) 

it follows that 

2 °° 

— / dk C n (fe) C m (fe) + / I C n (iK v ) CmjiKy) = S n ^ m . (21) 
7T J ,. 
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3 The inverse problem 

To take advantage of the algebraic analogue of the Marchenko method it is 
essential that there exist coefficients K n ^ m [independent of k] such that 

oo 

Cn(k)=J2K n)m f m (k). (22) 
m=n 

By analogy with Ref. ^T] assume that 

c n (k) = f n (k), n>N (23) 

[N specifies the order of a potential matrix]. If /jy and /jv+i are inserted 
[instead of respectively cn and cn+i] into Eq. (fTU|) for n = N, we obtain, in 
view of Eq. (jTJ, 

cjv-i(A;) = ( I^.jv-i/jv-iW + Pn iJV - a ^v] /jvW+ 

(24) 

.JV+l ~ /iV+l(&)] ) /&JV- 

Then, using the three-term recursion relation (fT3*|) with every n = N — 
1, . . . , 1 we obtain 

2JV-n-l 

C n (k)= Kn,mfm{k), 71 = 0, 1, . . . , N - 1 (25) 



[which in the limit N — > oo gives 

The coefficients i^ n>m are found from the completeness relation (|2ip. 
From the condition of the orthogonality of c n and every c m , m > n follows 
the condition of the orthogonality of c n and every f m , m > n, i. e. 



2 " 



- [ dkc n (k) f m (k) + Vc n (i^) / m (i«„) = 0, m > n. (26) 
o u 

Inserting the expansion of (|25jl in (f27)|) gives the system of linear equations 
in K n>m 

2N-n-l 

Kn, n Qn, m ^ t -^n, p Qp, 0, 771 !> 71. ("^) 

p=n+l 

Then, inserting Eq. (|2*5|l into (}2~T]) and putting n = m, we obtain, in view of 
Eq. (J26|) . the equation in K nyH 

-^n, n Qn, n ^ ] ^n,pQp,n I 1- (28) 

p=n+l / 
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Note that from Eq. ()27j) it follows that K n m , m > n are proportional to 
K n n . In Eqs. (J27j) and (J28|) Q n , m are defined from the scattering data by 



2 " 



). (29) 
o 

The elements a n and b n of the sought-for Hamiltonian matrix (jHJ) are 
related to K n ^ m by the equations 

rpl | n -|-l ^ j^n — 1, re Tn£ 

«n-J n> „-h Kn n J- n +l,n X„_i,„_i J n,n-1> 
w n — Kn _ ln _ 1 J n,n-1> ' ' ' 

ao is specified by the solutions Co and ci to Eq. |T0|) . 



4 A numerical realization 

To this point the assumption has been made that the phase shift 5g(k) is a 
continuous function of the wave number k that meets the conditions 



5*(oo) = 0, f k^Mk^dr < oo. 



In this case 5^ must satisfy stringent requirements. Indeed, the integrated 
function in r.h.s. of Eq. (j2~9"j) can be expressed in the form of a product g n g m 
of (real) functions 

g n {q) = S ne (q) cos5 e + C nl {q) sin^. (31) 

It is obviously that a sufficient condition to the convergence of the integrals in 
(j2HJ) is that functions g n are square-integrable. Notice that from (jl2)l follows 



i.e. the first term in (|3*Tj) decays exponentially at asymptotically large q. 
However, C n i grows exponentially with increasing q: 



C ne (q) ~ (-l)^ 1 ^-! ^ ^L q -K^^) e rr\ (33) 

q— +oo 7T 

This suggests that the phase shift 8g must decay rapidly enough to provide 
the convergence of the integral in r.h.s. of Eq. ()29|) . 



-(2n+e+2) q 2 /2 
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Actually the phase shift 5^ corresponding to the potential (JTJ of rank N 

m 

, s-n ^n-i,e(q) - ^iv(e) T N _i N S N t(q) 

tan<5 £ = -- — — t — --, 34 

as seen in Eqs. (I32J1 . (}3*3*j) . fulfills even more strict requirement 

rJV 7T (2A 4JV+21-3 -q 2 ,oc\ 

d < ^oo(iv-i)!r(iv + f + i) 9 ' 1 j 

Because of the restriction (}35j) on the phase shift 5f the potential (0) of finite 
rank iV generally is incapable to describe the scattering data on the infinite 
interval k € [0, oo). At most, we can set ourselves the task of constructing 
the potential (0) that describes the experimental phase shift 5e on some finite 
interval [0, ho], since generally 5g needs to be modified in the region k > k 
to provide at least the convergence of the integrals in Eq. (|29jh 

As an example we consider the s-wave scattering case. The "experimen- 
tal" phase shift 5i [dotted curve in figure 1] is that of the scattering on the 



potential given by straight well with the depth V : ^J^VqR — 1.5. The po- 
tential (JTJ) is sought for that describes the phase shift on the interval [0, k Q ] , 
k$R = 6 [in figure 1 crosses represent a modified phase shift]. The phase shift 
Sf corresponding to the resulting potential of rank = 6 and p = \R is 
shown in figure 1 [solid curve]. Notice that C n £ explodes exponentially with 
increasing q. Thus, the contribution from the region q > q = pk = 3 to the 
integral in Eq. (}2T?|) may become overwhelming [see figure 2 where go{q) 2 is 
plotted], with the result that the method fails. Matters can be improved by 
a transition to lesser p that shifts go to a region where C n e is not that large, 
or replacing 5e at q > q with a function that decays rapidly enough. 

In the second example, a scattering data on the potential with a bound 
state has been used as input. The phase shift Si [dotted curve in figure 3] 
corresponds to the s-wave scattering on a spherically symmetric potential in 



the form of straight well. The well parameter y |^ VqR = 2 determines the 
bound state with the energy E = —k 2 , kR = 0.638045 and the asymptotic 
normalization constant AdR 1 ^ 2 = 1.583324. k$ and p have been taken the 
same as in the first example [the modified phase shift is represented by crosses 
in figure 3]. It is well known that a phase shift does not depend on energy 
positions and asymptotic normalization constants of bound states. Thus, the 
inverse scattering problem in the presence of a bound state can be split into 
two steps. 

On the first step we focus on the describing the phase shift and, in spite of 
the whole of scattering data is used [see Eq. (|2H|) ]. do not seek to describe the 
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bound state with high degree of accuracy. The phase shift 5f corresponding 
to the potential ((TJ parameters, which together with kR and MR 1 / 2 are 
presented in the left half of Table, is shown in figure 3 [solid curve] . 

On the second step, to improve the description of the bound states we 
use the relationship (fTTIJ) between the poles and residues of the S-matrix 
and the characteristics of the bound states [see e.g. Ref. [S]]. Here, the 
smallest eigenvalue A and the corresponding eigenvector component Zjv-i,o 
associated with the bound state are found from the system 

N-l 

E Z n-u = 1, (36) 

3=0 

P n (-k p /2) = — £ , (37) 

N—l, N C N ' e {lKp) 

Cn-ij^P) ~ Pw(-KV/2)T;_ liJV 4;j(kp)_ = P M 2 (38) 



d {C^! h M ~V N (q 2 /2)Tf,_ ltN C^(q)} 



dg 



q=\Kp 



Eqs. f!37jl . (jHHjl are derived from the S-matrix formula for the potential ((TJ) 

13 



s ,,\ = C { N l i e (q) - Vn(£) t n-i,n Cn,\{q) -, f) 
C^e(q)-V N (e)Tt N CtiUq) 



and Eq. (JTJJ). Notice that the component Zjv-i,jv-i corresponding to the 
leading eigenvalue Xn-i is involved to meet the normalization condition (jHEJ)- 
The phase shift is scarcely affected by changing the parameters {Ao, Zjv-i,0j 
Zn-i,n-i} from the initial values obtained on the first step to the ones that 
are evaluated from Eqs.(|nEJ) - J3EJ)- The potential parameters, which provide 
the correct values of kR and A4R 1 / 2 , are presented in the right half of Table. 



5 The Laguerre basis 
5.1 Preliminaries 

For simplicity's sake we restrict our consideration to the scattering of neutral 
particles. However, the resulting equations still stand in the presence of 
the repulsive Coulomb interaction. The potential sought is given by the 
expression 

*2 N-l _ _ 

Vi(r, r') = — J2 ^ n (x)V n ,J m (x') (40) 

Z P n, m=0 
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where the functions <j) n 

*M = r(B + M + i)i *M (41) 

are bi-orthogonal to the base Laguerre functions <p e n : 

cj) e n (x) = (2br) i+1 e- br L 2 n e+1 (2br), (42) 

i.e. 

drfyx) <fi e m (x) = 6 njm . (43) 

Here, b is the scale parameter: x = br. 
The coefficients u n of the expansion 

oo 

1>(k, r) = Y,Mk)<p £ n (x) (44) 

n=0 

of the Schrodinger equation regular solution ip(k, r) satisfy the system of 
equations 

(h° n ,m + V n,m)um(k) = k 2 A e nm U m (k), 71 = 0,1,.... (45) 

Here, \\h^ m \\ is the symmetric tridiagonal matrix of the reference Hamilto- 
nian © calculated in the basis (02)) [Q: 



<n = ^(n + Hl), 
ft - ft - fe (n+2£+2)! n - 1 

"'n,n+l ~~ "'n+l.n ~ u 2n\ ' ' ' 



(46) 



|A^ m || signifies the basis-overlap matrix 

oo 

Ai m = ( dr<\> l n {x)<\> l m {x) (47) 



which is also of Jacobi form: 



A* - A , - n - 1 



(4f 



The asymptotic behaviour of the coefficients u n (k), k > 0, as n — > oo is 
given by the following expression: 

u n (k) = w n (k) = i - S(*)C$(*)] (49) 
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where the functions [7j 

n! 

(n+e+l)\ (2sinC) 



obey the inhomogeneous "free" equation 

J n,m(k) C£? e (k) = 5 n>0 , n = 0, 1, .... (51) 

Here, ||J^ m (fc)|| = m — k 2 A e n ^ m \\ is the so-called J- matrix. S n j are the 
solutions of the system of equations 

Jn, m (k)S mJ (k) = 0, n = 0, 1, (52) 
/' C2 sinH m 

^n,/(fc)= 2(2£+1) , (-O w 2^i(~n^ + l; 2^ + 2; 1-T 2 )- (53) 
It can easily be shown that the completeness relation for the functions S n ^ 



of Ref. [in] can be rewritten as 

2 °° 

- / dkS nJ (k)A e n , m S mJ (k) = 5 n , n , (54) 

7T J 


The coefficients M n (i/t J/ ) of the expansion of the bound state normalized 
wave function ip u (r) with the energy — k 2 have the following asymptotic be- 
haviour 

u n (iK v ) = w n (iK u ) = M v i l C^j(iK v ) (55) 

asn-> oo. 

oo 

Notice that the sine-like J-matrix solutions S(r) = J2 S n j(k) 4>n{x) and 

n=0 

the cosine-like one C(r) = E C e(k) <Pn(x), where C n Ak) = I (C„ l{k) + 

C„~j(fc)V have the asymptotic behaviour (fTHjl. 

The the completeness relation (|2T?|) is transformed into 

2 r 

- / u n (k) A , m u rn (k) + u n (iK u ) A n , u m {w v ) = 5 nn >. (56) 

71 J 
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5.2 Inverse problem 



In the framework of the J-matrix version^] of the inverse scattering problem 
the spectral parameter set {Xj, Z^-i ,j}j= is obtained from the scattering 
data of the truncated Hamiltonian matrix of order N in the orthogonal basis 

2— it 



Dt 



d n , n > m, 



n < m, 



2bn\ 



(n + 2£ + 2)\ 



(57) 



i. e. 



d e n (2brY +1 e- br L 2 n e+2 (2br). 



Clearly the set {Xj, Zjv-i, j}^^ determines a tridiagonal Hamiltonian matrix 

of order N in any orthogonal basis x n = Em=o P% m where \\P* m \\ is an 
arbitrary orthogonal (N x iV)-matrix of the form 



(58) 



ipe 



0,0 



P 



0,N-2 



P 



N-2,0 



P 



V o 



N-2,N-2 





o \ 



1 ) 



(59) 



Let us assume that \\P~ m \\ is the orthogonal transformation matrix that 

f p 1 N-l r f -i N-l 

performs the change from jv^nj, to the new basis |Xnj m "which the 
kinetic energy operator truncated matrix is tridiagonal. To perfect the anal- 
ogy with the oscillator basis case, let us denote the kinetic energy operator 

2u - - . . . - . r ^N-l 

7? 



H° ((7j) tridiagonal matrix in the basis {xl} Q by ||T^ m ||. The sought for 
Hamiltonian j^H matrix ||/i nim || of order N is presumed to be of a Jacobi 

r . > N-l 

form (p) in the basis jXnj • 

Thus, the first iV — 1 of the wave function ip(k, r) expansion coefficients 
in the combined basis set {{x n }n=oi {01}^tv} obey the equations 



a c (k) + b x cx{k) = k 2 c (k) 

On C n —\ 

(k) + 

d n c n (k) + b n+ i c n+ i(k) = k 2 c n {k), n 



1, ... N -2. 



(60) 



It is easy to check that a sufficient condition that the algebraic version of 
the Marchenko method be applicable for the construction of the tridiagonal 
Hamiltonian matrix (R3) is that 



T 



ra+l, ra+1) 



+i 



T 



n, n+1 ' 



for n 



M,...,N 



2, M 



rfi- (si) 
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If N is odd, to (J5TJ) must be added the constraint that cim = M . In this 
case cm+i = /m+i> c m = Im, where f n satisfy the "free" equations 

Ti, n -J n -i(k) + T l n J n (k) + Tl n+1 f n+1 (k) = k 2 f n (k), n = 1, . . . , N - 2, 

(62) 

and we obtain for n < M — 1 

N-n-l 

C n {k) = K n,mfm{k). (63) 

m=n 

Notice that in going from the initial Laguerre basis {^i} to the com- 
bined basis set {{x£}n=o\ Wn)n=N] the submatrices ||fr° im || J~L and 
ll^lmllnm=o are transformed into ||7^ m ||^~l an d the identity matrix of or- 
der N respectively In addition, the elements h%_ 1N , A e N _ 1N and h% N _±, 
A e N n-1 are multiplied by d e N _ v The rest of the (infinite) matrices ||/i° m || 
and H-A* m || is unaltered. It thus follows that f n must satisfy (in addition to 
(|52)l) the equations 

T N-i,N-2fN~2(k) +T%_ ltN _ 1 f N -i(k) +d e N _ 1 f N _ 1 N f N (k) = k 2 f N „x(k), 

(64) 

d^-iJff,N-i{k)f N -i{k) + Jff,N( k )M k ) + J N,N + i(k)fN + i(k) = 0. (65) 

Jn,m(k)fm(k) = 0, Tl = N + 1, . . . . (66) 

Putting /„ = S n ^ for n > N and /jy-i = Sn~i,i/ d l N _ x [in view of the 
equation (JBTIj) and (|63j) . respectively], from Eq. by using the tree-term 
recursion relation (|62jl we obtain the coefficients S n / with n = 0, . . . ,N — 2. 
Similarly, setting f n = C^} for n > N and /jv-i = ^N-iil ^N-x-t we obtain 
the coefficients with n = 0, . . . ,N — 2. From the Wronskian-like relation 
(see e.g. [TH] ) 

- CgW^vC*)) = fc, n > (67) 

it follows that 

It-M (^tl#)«M*0 - C$(k)S n+1>e (k)) =k,0<n<N-2. (68) 

Besides, since the system of equations ()52j) in <S n ^ is homogeneous, the sets 
{^n/j^Lo anc ^ {"-w} are connected by a linear transformation and there- 
fore S n< £ also satisfy the homogeneous equation 

T o VV(*0 + T^St^k) = k 2 S 0/ (k), (69) 
13 



whereas obey the inhomogeneous one 

T^{k) + T* tl c[f(k) = ec$(k) + (70) 

Thus, the two sets, {<S n ,e} n=0 and {C n ,e} n=(j , C n ,e = + C$), are 

"free" independent respectively sine-like [S n ,i = <Sn,e, n> N] and cosine-like 
[C& = C$, n>N] solutions to Eqs. ©-fH (see e.g. [7j). 

From the above discussion it follows that to obtain f n with < n < N—l, 
which are involved in Eq. (JBHJ), we can place /jy = w^, /n-i = Wn-i/^n-xi 
where w n are defined by (f4T)J) . Then, inserting this /at and Jn-x in Eq. (fMJ) 
gives /n-2- Once /at_2 and /jv-i are known, /„ for n = iV — 3, . . . , are 
obtained by using the tree-term recursion relation ()62|). K n ^ m are determined 
by the equations (j2Zj) and (j^Hj) [in which the upper limit in the sums is equal 
to N — n — 1]. The expressions for {a n }, {&„} are the same as (J3TJ|) . 

Qn,n' with n < iV — 1 in Eqs. (f2"Tj) and (f2"KJ) are defined [in view of the 
overlap matrix form in the combined basis] by 

Qn,n> = - I fn{k)jAF) +J2fn(^u)jA^uj, Tl' < N - 2, (71) 
7T J 



Qn,N-i = ~ J fn(k) [wN-iify/db^ + d l N _ 1 A e N _ ltN w N {k)} + 



+ E/n(i« 1 /) [wjV-l(i«*0/<%-i + d N-l^N-l,iV w Jv(iK«/)], 

(72) 

oo 

2 r 

Qn,n' = - / f n (k) A e n , m w m (k) + f n (xK v ) A n , m w m (xK v ) , n' > N. (73) 

7T J ,, 



Notice that at large k, as seen in Eq. O, C^{k) ~ jfe*. Thus, as in the 
case of the oscillator basis, we should restrict ourselves to the description of 
the scattering data on a finite energy interval, beyond the boundary of which 
the phase shift needs generally to be modified to provide the convergence of 
the integrals in Eqs. (|7T|) - (f73)l . 



6 Conclusion 

In the potential scattering case the finite-difference approach and J-matrix 
method share the tridiagonal representation of the Hamiltonian. The analogy 
can be carried over to the inverse scattering problem formalism. Here, the 
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J-matrix version of the Marchenko equation algebraic analogue is formulated 
and its numerical realization features are considered. The merit of JME is 
that it is free from a parameter fit inherent in the previous J-matrix inverse 
scattering approach El HI E] • We also construct a tridiagonal Hamiltinian 
matrix of some order M in an orthogonalized Laguerre basis; in doing so 
it is sufficient to tridiagonalize the matrix representation of the reference 
Hamiltonian H° in the finite orthogonal basis subset of size N = 2M. As 
has been shown in Ref. [Hj, in the two coupled-channel case without threshold 
the sought-for interaction matrix may be of a "quasi-tridiagonal" form. On 
this assumption JME can be easily expanded to multichannel scattering. 
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Table 







N = 7, p 


- R 






kR = .6512647458, 


kR = 


6380449999, 




MR 1 ' 2 = 1.6017576599 


MR 1 / 2 = 


= 1.5833238674 


3 




\ 3 


ZN-l,j 


A, 





0.0356514517 


-0.0381260178 


0.0362075259 


-0.0353279575 


1 


0.1482712147 


0.4605384282 


0.1482712147 


0.4605384282 


2 


0.2309801539 


1.3246452781 


0.2309801539 


1.3246452781 


3 


0.3094585382 


2.5044702689 


0.3094585382 


2.5044702689 


4 


0.4084394267 


4.1865934000 


0.4084394267 


4.1865934000 


5 


0.5275945630 


6.7063360348 


0.5275945630 


6.7063360348 


6 


0.6184249465 


10.1425219887 


0.6183926387 


10.1425219887 
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